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1  Introduction 


There  are  two  main  reasons  for  transverse  gas  injection  into  supersonic  flows.  One  is  to  mix 
fuel  into  a  supersonic  air  stream  for  supersonic-combustion  ram  jet  propulsion,  where  the 
emphasis  is  on  quality  of  mixing.  The  other  is  to  cause  a  reaction  force  on  a  vehicle  traveling 
at  supersonic  speed,  or  in  a  supersonic  nozzle  for  thrust-vector  control,  where  the  emphasis 
is  on  effectiveness  and  undetectability  of  vehicle  control.  In  both  cases  good  design  depends 
crucially  on  the  capability  to  compute,  and  thus  to  predict,  the  associated  complex  turbulent 
flow  field. 

At  the  outset  of  the  program,  the  GALCIT  supersonic  Ludwieg  Tube  (LT)  was  within  six 
months  of  being  completed.  This  facility  provides  a  flow  at  Mach  number  2.3  in  a  test  section 
of  200x200  mm  cross  section  for  a  duration  of  80  ms.  It  is  very  inexpensive  to  operate.  A 
secondary  Ludwieg  Tube  for  synchronized  gas  injection  into  the  supersonic  flow  was  under 
development. 

At  the  same  time,  the  development  of  digital  defocusing  particle  imaging  velocimetry 
(DDPIV)  by  the  Gharib  group  at  GALCIT  for  use  in  low  speed  water  flows  had  made 
it  possible  to  obtain  three-dimensional,  three- component,  instantaneous  measurements  of  a 
complete  velocity  field.  The  combination  a  suitable  variant  of  this  tool  for  supersonic  flow 
with  the  LT  promised  to  make  the  acquisition  of  such  flow  fields  highly  productive  and  inex¬ 
pensive.  Comparing  the  method  with  regular  or  even  Stereo  DPIV,  where  the  data  from  one 
run  would  be  a  velocity  vector  field  in  just  one  plane,  so  that  many  runs  would  be  required 
to  get  a  velocity  field  on  a  three-dimensional  grid  (and  which  therefore  would  not  be  tem¬ 
porally  correlated),  the  productivity  is  increased  by  a  factor  of  20  or  so,  and  the  value  of  the 
temporal  correlation  of  the  vector  field  is  huge.  With  an  expected  turnaround  frequency  of 

the  LT  of  6  runs  per  day,  and  at  an  operating  cost  of  $  400  per  day,  this  was  a  very  attractive 
proposition. 

More  sophisticated  turbulence  computations  are  being  made  accessible  by  rapid  develop¬ 
ments  in  computational  power  and  data  storage  of  modern  machines  every  year.  More 
sophisticated  turbulence  models  are  being  developed  at  the  same  time.  Computational  fluid 
dynamicists  are  asking  for  more  complete  and  more  highly  resolved  data  than  are  currently 
available  for  transverse  injection  into  supersonic  flow. 

The  motivations  for  the  work  was  thus  to  investigate  the  flow  field  of  transverse  injection 
into  a  supersonic  flow  with  fine  spatial  resolution  of  the  instantaneous,  three-dimensional, 
three-component  velocity  field,  the  productivity  of  the  method  permitting  coverage  of  large 
ranges  of  the  parameter  space. 

In  order  to  be  able  to  develop  the  instrument,  a  separate  proposal  for  a  DURIP  Grant 
was  made.  This  was  successful,  and  the  instrument  was  built  in  its  first  form  under  the 
DURIP  Grant.  It  rapidly  became  clear  during  further  development  and  calibration  of  the 
instrument,  that  the  step  from  low  speed  water  flows  to  supersonic  flows  at  Mach  2.3  in  air 
involves  major  development.  The  original  statement  of  work  (see  Appendix  A)  thus  turned 
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out  to  be  much  too  ambitious.  In  the  three- year  program  we  can  therefore  only  claim  to 
have  achieved  the  development  of  the  instrument  to  the  stage  where  it  can  be  applied  with 
success  in  supersonic  air  flows.  In  order  to  achieve  this,  the  instrument  had  to  be  significantly 
modified,  and  entirely  new  software  had  to  be  developed.  The  major  part  of  this  final  report 
describes  this  work. 

Before  getting  to  that  part,  the  next  section  describes  the  Ludwieg  Tube  and  the  secondary 
Ludwieg  Tube. 


2  The  supersonic  Ludwieg  Tube 

A  Ludwieg  Tube  consists  of  a  long  tube,  filled  initially  with  high-pressure  gas,  that  is 
connected  via  a  nozzle  and  test  section  to  an  evacuated  dump  tank.  The  dump  tank  and 
test  section  are  separated  by  a  diaphragm  which  is  ruptured  to  initiate  the  flow.  During  the 
time  it  takes  for  the  expansion  wave  generated  by  the  diaphragm  rupture  to  propagate  to 
the  far  end  of  the  Ludwieg  tube  and  back,  the  flow  into  the  nozzle  is  perfectly  steady  and 
uniform  except  for  the  boundary  layer  on  the  tube  wall,  which  grows  with  time  and  becomes 
turbulent  after  a  certain  time. 

A  donation  from  the  John  and  Anna  Wild  estate  enabled  us  to  design  and  build  a  Ludwieg 
Tube  to  replace  the  small  continuous-flow  supersonic  tunnel  that  was  scrapped  as  part  of 
the  Guggenheim  building  rehabilitation.  The  LT  has  been  operational  since  2001.  The  total 
cost  of  the  facility  is  $  500,000.  The  first  nozzle  of  the  LT  is  designed  for  a  Mach  number 
of  2.3  with  a  test  section  cross  section  of  200x200  mm.  The  minimum  free-stream  viscous 
length  scale  v^/U^  =  10  nm,  giving  a  maximum  Reynolds  number  on  a  length  of  200  mm  of 
20x10®.  The  tube  has  a  diameter  of  300  mm  and  a  length  of  17  m  honed  to  a  smooth  finish. 
This  provides  a  theoretical  test  duration  of  100  ms.  The  starting  process  takes  approximately 
20  ms,  so  that  the  clean  test  duration  is  80  ms,  or  200  flow  times  for  a  200  mnn  model.  The 
short  test  time  is  the  reason  for  the  inexpensive  operation  of  the  LT.  A  perspective  view  of 
the  facility  is  shown  in  Fig.  1. 

A  particular  feature  of  the  GALCIT  LT  is  the  method  by  which  the  tube  boundary  layer 
is  removed.  From  experimental  results  obtained  by  Ghun  (private  communication)  it  has 
become  clear  that  the  growth  of  the  boundary  layer  in  a  LT  and  the  transition  to  turbulent 
flow  in  this  boundary  layer  can  provide  serious  impairment  of  the  flow  quality,  and  this 
makes  it  unattractive  to  extend  the  tube  beyond  approximately  40  diameters.  Even  then, 
it  is  desirable  to  suck  off  the  tube-wall  boundary  layer.  This  is  done  by  providing  a  flow 
path,  parallel  to  that  through  the  nozzle,  through  an  annular  throat  just  upstream  of  the 
contraction,  through  a  manifold  and  four  tubes  to  the  space  upstream  of  the  diaphragm,  see 
Fig.  2.  So,  when  the  diaphragm  is  ruptured,  an  expansion  wave  passes  both  through  the  test 
section  and  nozzle  into  the  tube,  and  through  the  tubes,  manifold  and  annular  throat  into 
the  tube.  This  sets  up  both  flows.  As  the  tube-wall  boundary  layer  grows,  the  rounded  nose 
on  the  inside  of  the  annular  throat  permits  the  stagnation  point  to  move  without  causing 
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Figure  1;  Perspective  view  of  Ludwieg  tube,  nozzle,  test  section  and  dump  tank 


separation.  A  computation  of  this  flow  indicates  that  this  method  of  boundary-layer  removal 
starts  the  flow  successfully.  The  region  around  the  annular  throat  is  instrumented  with 
pressure  transducers  in  order  to  monitor  the  performance  of  the  removal  system.  During 
the  calibration  of  the  LT,  these  transducers  showed  that  the  flow  downstream  of  the  annular 
throat  is  supersonic,  thus  ensuring  that  no  noise  can  enter  the  main  flow  from  the  suction 
manifold. 


2.1  LT  nozzle  calibration 

The  two-dimensional  Mach  2.3  nozzle  was  designed  by  J.  J.  Korte  of  NASA  Langley  RC. 
The  complete  nozzle  flow  was  calibrated  [1]  by  using  the  Mach-wave  method  in  which  weak 
disturbances  are  generated  at  opposite  nozzle  walls  by  sticking  strips  of  thin  adhesive  tape 
to  the  walls.  The  weak  Mach  waves  generated  by  the  tapes  can  be  visualized  by  the  schlieren 
technique.  The  angles  of  the  intersecting  Mach  waves  give  the  local  Mach  number  and  flow 
direction.  A  composite  view  of  the  schlieren  images  is  shown  in  Fig.  3,  and  the  resulting 
Mach  number  distribution  is  compared  with  the  design  curve  in  Fig.  4 
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Figure  2:  Longitudinal  section  through  nozzle  and  test  section.  To  the  right,  this  is  joined  by  the 
17m  long  300  mm  diameter  tube,  and  to  the  left  by  the  diaphragm  rupture  system  and  dump  tank. 
The  light  blue  and  yellow  parts  at  the  right  form  the  annular  throat  of  the  boundary-layer  removal 
system.  The  magenta  part  shows  the  transition  from  the  circular  section  to  the  rectangular  throat. 
The  four  tubes  connecting  the  two  green  manifolds  of  the  boundary-layer  removal  system  are  not 
shown. 


Figure  3:  Composite  schlieren  picture  of  the  whole  nozzle  flow,  made  up  of  seven  separate  images 


2.2  The  injection  system 


In  order  to  be  able  to  control  the  state  of  the  boundary  layer  on  the  surface  from  which 
the  gas  injection  takes  place,  a  flat  plate  is  constructed  that  is  located  parallel  to  and  a 
small  distance  below  the  test  section  roof.  This  enables  the  tunnel-wall  boundary  layer  to 
be  diverted.  The  injection  jet  is  generated  with  a  separate  small  injection  Ludwieg  Tube 
(ILT)  of  50  mm  internal  diameter  and  10  m  length.  This  feeds  the  injection  nozzle  via  a 


6 


Figure  4:  Mach  number  distribution  obtained  from  measurements  of  the  intersection  angles  of 
the  waves  in  Fig  3.  The  diamonds  represent  measurements  of  the  intersection  of  the  weak  shocks. 
The  squares  represent  measurements  from  intersections  of  the  leading  characteristics  of  the  trailing 
expansion  waves.  The  full  line  is  the  centerline  Mach  number  distribution  obtained  by  J.  J.  Korte 
in  his  design  of  this  nozzle,  which  takes  account  of  the  boundary  layer  displacement  thickness. 


mylar  diaphragm  that  is  ruptured  by  discharging  a  capacitor  through  a  tungsten  wire  that 
is  clamped  onto  the  mylar  diaphragm.  The  local  heating  of  the  mylar  diaphragm  causes  a 
disk-like  piece  of  the  diaphragm  to  be  cut  out.  This  permits  accurate  synchronization  of  the 
injection  flow  with  the  main  flow.  The  maximum  injection  nozzle  exit  diameter  is  12  mm. 
This  large  diameter  makes  the  size  of  the  flow  field  of  interest  in  the  longitudinal  vertical 
plane  approximately  150x200  mm.  The  maximum  jet  Reynolds  number  is  4x10®. 

The  injection  is  timed  to  start  immediately  after  the  starting  process  of  the  main  nozzle  flow 
is  over  and  the  boundary-layer  removal  system  has  reached  steady  state.  This  will  be  no  later 
than  20  ms  after  main  diaphragm  rupture.  The  flow  duration  of  the  ILT  is  approximately 
60  ms. 

A  schlieren  image  of  the  flow  with  the  ILT  in  synchronized  operation  with  the  LT  is  shown 
in  Fig.  5. 
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Figure  5;  Schlieren  image  of  the  jet  in  cross  flow.  The  exposure  time  is  too  long  to  resolve  the 
turbulent  structure  in  this  image 

3  The  three-dimensional  particle  tracking  velocimetry 
system 

3.1  Introduction 

During  the  last  few  decades,  flow  experiments  where  three-dimensional  features  are  of  interest 
have  become  more  common.  This  evolution  brings  the  need  for  a  full  three-component,  three- 
dimensional  velocity  field  measurement  technique..  For  that  purpose,  some  researchers  have 
developed  methods  using  a  single  camera  [2]  or  several  cameras  (generally  3  or  4  using 
the  epipolar  technique)  [3,  4,  5].  Of  particular  interest  is  the  work  done  by  F.  Pereira  et 
al.  [4,  6]  in  developing  a  digital  defocusing  particle  image  velocimetry  (DDPIV)  system. 
Their  system  uses  three  individual  CCDs  placed  at  the  vertices  of  an  equilateral  triangle, 
with  all  three  cameras  having  overlapping  fields  of  view.  The  test  section  can  be  placed 
anywhere  in  this  overlap  region.  Every  particle  in  the  test  section  will  be  projected  onto 
each  CCD.  If  the  CCD  images  are  overlaid,  the  three  images  of  a  particle  will  be  located 
at  the  three  corners  of  a  near-equilateral  triangle.  By  performing  triangulation,  the  three- 
dimensional  position  of  a  particle  can  be  calculated.  Sub-pixel  accuracy  of  a  particle  image 
location  is  obtained  by  having  the  images  unfocused,  which  drastically  improves  the  overall 


resolution  of  the  system.  This  design  was  patented  by  Gharib  et  al,  and  is  assigned  to 
the  California  Institute  of  Technology  (U.S.  Patent  #6,278,947).  Currently,  this  technology 
is  being  developed  for  commercial  use  by  VioSense  Corporation  of  Pasadena,  California. 
F.  Pereira  et  al.  also  developed  a  sophisticated  software  toolkit  which  allows  users  to 
analyze  the  raw  CCD  images  and  transform  them  into  a  three-dimensional  velocity  field. 
The  DDPIV  system  has  been  applied  to  bubbly  flows. 

The  current  authors  first  revisited  some  of  the  design  choices  made  by  F.  Pereira  et  al,  and 
improved  upon  some  of  them.  Also,  since  the  source-code  of  F.  Pereira  et  al  is  proprietary, 
a  completely  new  software  package  was  developed  and  is  significantly  different  from  the 
DDPIV  software  package.  The  current  authors  will  investigate  supersonic  air  flows,  whereas 
F.  Pereira  et  al  investigated  low-speed  water  flows.  Of  course,  the  software  of  both  F.  Pereira 
et  al  and  the  current  authors  share  some  of  the  same  theoretical  foundation;  however,  much 
of  the  new  software  is  significantly  different  and  substantially  improved.  These  improvements 
provide  increased  accuracy  of  results.  F.  Pereira  et  al  chose  to  use  a  statistical  calculation  for 
the  velocity  field;  whereas,  the  current  authors  pursue  direct  particle  tracking.  The  system 
developed  by  the  current  authors  is  referred  to  as  three  dimensional  particle  triangulation 
velocimetry  (3DPTV). 


3.2  Principle 

Both  the  hardware  setup  of  F.  Pereira  et  al  for  DDPIV  and  the  hardware  setup  of  3DPTV 
consist  of  a  camera  comprising  three  separate  CCDs,  which  all  view  the  same  test  volume 
from  different  angles.  The  front  plate  of  the  camera  consists  of  a  mask,  three  lenses,  and 
three  apertures,  behind  each  of  which  there  is  a  CCD.  The  need  for  more  than  two  CCDs 
and  suggestion  that  with  three  cameras  an  equilateral  triangle  is  the  optimum  configuration 
is  discussed  by  Maas  et  al  [3].  In  the  DDPIV  setup,  the  center  of  a  target  plate  is  focused 
onto  the  center  of  each  CCD.  Figure  6  shows  this  for  one  of  the  CCDs.  In  this  arrangement, 
the  image  of  the  test  volume  does  not  fill  the  frame  of  the  CCD.  In  order  to  make  full  use 
of  all  the  pixels  of  the  CCD,  3DPTV  uses  offsets,  in  which  the  CCDs  are  moved  away  from 
the  lens  and  away  from  the  camera  axis  as  shown  in  Figure  6.  Of  course,  full  use  of  the 
CCDs  could  also  be  achieved  by  placing  the  target  plate  in  the  test  volume,  but  this  has 
some  disadvantages,  as  is  discussed  later.  The  image  of  each  point  on  the  target  plate  will 
appear  at  the  same  place  on  each  of  the  three  CCDs,  once  the  correction  for  the  offsets  has 
been  made. 

In  practice  it  is  impossible  to  place  the  CCDs  in  their  exact  theoretical  locations;  therefore, 
careful  alignment  followed  by  calibration  is  required.  In  order  to  do  this,  a  well  designed 
target  is  needed  in  addition  to  calibration  software.  For  this  purpose  the  authors  designed  a 
target  with  100  equally  spaced  dots  per  square  inch.  With  this  target  placed  at  the  target 
plane,  the  CCDs  should  be  placed  such  as  to  focus  and  match  each  of  the  three  images  as  well 
as  possible.  The  calibration  software  compares  the  measured  and  known  dot  locations,  and 
then  creates  a  mapping  to  correct  for  errors  in  CCD  placement  as  well  as  lens  aberrations 
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Figure  6:  Maximum  CCD  utilization  using  offsets. 
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Once  the  camera  has  been  aligned  and  calibrated,  the  images  of  a  particle  that  does  not  lie 
on  the  target  plane  will  appear  at  the  vertices  of  an  equilateral  triangle  on  the  overlaid  CCD 
images.  The  formation  of  triangles  on  the  overlaid  images  is  shown  in  Figure  7.  Clearly,  it 
is  only  a  matter  of  ray  tracing  to  transform  a  triangle  on  the  CCDs  to  a  particle  location. 

Particle  tracking  is  conceptually  a  simple  task.  Once  all  of  the  three  dimensional  positions 
are  known,  all  that  is  required  to  determine  the  three  components  of  the  velocity  of  the 
particles  is  to  find  corresponding  particles  in  two  successive  frames  separated  by  a  known 
short  time  interval. 


3.2.1  3DPTV  Design 

The  design  of  the  3DPTV  system  involves  many  related  parameters.  Choosing  the  best 
parameters  for  a  set  of  particular  requirements  is  truly  a  design  exercise.  The  governing 
equations  used  in  this  design  are  now  presented.  A  key  assumption  in  the  design  is  that 
the  medium  surrounding  the  camera  has  a  constant  refractive  index.  The  fact  that  this 
assumption  is  not  true  means  that  the  system  will  be  used  in  an  off-design  regime;  however, 
for  the  current  parameters  this  leads  to  only  a  very  small  deviation  from  the  optimum. 

Blur  Radius  Estimation 

The  blur  radius  has  a  big  influence  on  the  results  one  can  obtain  with  both  DDPIV  and 
3DPTV.  A  sufficiently  large  blur  radius  is  required  to  get  sub-pixel  accuracy;  however,  the 
larger  the  blur  radius,  the  more  often  particle  images  will  overlap,  thereby  making  it  more 
difficult  to  process  the  data.  An  estimate  of  the  blur  radius  can  be  calculated  by  considering 
the  quasi-one-dimensional  problem  of  a  single  lens  of  focal  length  /  and  a  single  CCD  as 
shown  in  Figure  8.  The  distance  between  the  lens  and  CCD,  in  the  2-direction,  is  dc  and 
the  lens  has  an  aperture  of  radius  r^.  The  setup  was  constructed  so  that  a  particle  placed  at 
a  distance  L  from  the  camera  will  be  focused  on  the  CCD.  If  a  particle  is  at  a  distance  2,  it 
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Test  Volume 


Figure  7:  Projection  of  a  particle  onto  the  CCDs,  with  appropriate  offsets  taken  into  account. 


Figure  8:  Quasi  1-D  blur  radius  evaluation. 
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will  be  focused  away  from  the  CCD,  at  a  distance  fz/{z  -  /).  Geometric  relations  lead  to 
an  estimation  of  the  blur  radius,  r^, 


/  ^ 

The  rate  of  change  of  the  blur  radius  to  a  movement  in  ^  is 


drh  _  fL 

dz 

The  data  processing  program  has  to  find  the  centroids  of  all  the  Gaussian  intensity  distribu¬ 
tions  of  the  blurred  particle  images.  The  program  does  not  know  the  blur  radius  of  a  particle 
image  a  priori  since  the  particle  location  is  unknown.  The  performance  of  the  program  is 
improved  by  operating  with  a  small  range  of  blur  radii.  The  value  of  the  average  blur  radius 
and  range  of  blur  radii  are  parameters  that  must  be  considered  when  choosing  both  L  and 
/.  One  should  also  note  that  a  nonzero  blur  radius  is  required  to  achieve  sub- pixel  accuracy; 
therefore,  the  focal  plane  can  not  be  located  within  the  test  volume. 


Alignment 


As  explained  before,  the  camera  is  composed  of  three  GGDs  placed  at  the  vertices  of  an 
equilateral  triangle.  Let  the  coordinates  of  the  center  of  the  GGD  be  P?.  Each  of  the 
GGDs  views  the  test  volume  through  its  own  lens,  of  focal  length  /,  which  is  located  at  Pl. 
An  object  placed  at  the  focal  distance  of  the  camera,  L,  will  look  exactly  the  same  on  each 
GGD,  after  appropriate  offsets  have  been  taken  into  account,  and  will  be  focused.  Therefore, 
the  distance  in  the  ^-direction,  between  the  GGD  and  the  optical  center  of  the  lens,  dc,  is 
given  by 


dc  = 


fL 


(1) 


where  r  is  the  radius  of  the  circle  on  which  the  center  of  the  lenses  are  located.  One  possible 
lens  configuration  is 


Pl  =  (0,r,0) 


Finally,  if  the  camera  is  perfectly  aligned,  dhe  positions  of  the  GGD  centers  for  this  lens 
configuration  will  be 
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The  purpose  of  alignment  and  calibration  is  to  compensate  for  the  placement  of  the  CCDs 
as  to  agree  with  the  above  equations.  This  is  a  critical  step  to  ensure  the  accuracy  of  the 
measurements. 

SD  Coordinates  and  Sensitivity 


Suppose  that  a  particle  is  at  position  {x,  y,  z).  Using  geometric  optics  it  is  possible  to  find 
the  pixel  coordinates,  U,  of  the  particle  image  on  CCD  i: 


where 


(-!). 


and  where  S  is  the  pixel  size.  Using  these  coordinates,  the  side  length,  s,  of  the  equilateral 
triangle  formed  by  overlaying  the  three  CCD  images  is  given  by 


(2) 


It  is  now  easy  to  obtain  the  sensitivities  of  a  3DPTV  camera: 


do, 

~d^~~f  ~Yz 

ds  \/3r  , 

Tz  ~  ('i) 

The  sensitivity  in  the  x-  and  y-directions  are  the  same;  however,  in  most  cases,  the  sensitivity 
in  is  much  lower,  about  10  times.  It  is  interesting  to  note  that  the  ratio  of  these  sensitivities 
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depends  only  on  r  and  z, 


dli.z 

dx  _  ^ 

5z 

Equations  1,  3  and  4  show  that  higher  /  leads  to  better  accuracy.  Of  course,  /  is  constrained 
by  space  limitations  and  by  the  need  for  sufficient  depth  of  field.  Increasing  r  increases  the 
accuracy  in  2:  and  decreases  the  accuracy  in  the  other  directions.  Decreasing  the  distance 
between  the  camera  and  the  test  volume  increases  all  the  accuracies.  However,  a  refractive 
index  change  along  the  light  path,  caused  for  example  by  the  presence  of  a  glass  window 
or  non-uniform  surrounding  medium,  introduces  a  difficulty  if  r  is  not  small  compared  to  z. 
This  difficulty  is  discussed  further  in  the  next  paragraph.  Finally,  in  most  cases,  decreasing 
L  increases  the  accuracy.  However,  the  blur  radius  range  limitation  does  not  allow  us  to 
reach  the  highest  sensitivities. 


3.2.2  Varying  Media  Influence 

In  most  of  the  results  presented  here,  a  glass  window  has  to  be  taken  into  account.  Although 
the  sensitivities  are  not  significantly  changed,  the  triangles  are  no  longer  equilateral.  This 
poses  an  additional  problem  to’  the  data  processing  program.  To  keep  the  triangles  nearly 
equilateral,  the  light  path  from  a  particle  to  the  CCD  has  to  have  a  small  angle  of  incidence 
at  the  glass.  This  leads  to  two  important  conditions.  The  first  is  that  the  camera  has  to 
look  perpendicularly  through  the  glass  and  the  second  is  that  the  ratio  rjz  has  to  remain 
much  less  than  1. 


3.2.3  Current  Camera 

The  current  camera  is  designed  to  investigate  a  50  x  50  x  50  mm^  test  volume.  The  CCDs  are 
UNIQ  UP-1830-10  cameras.  Each  sensor  is  composed  of  1024  x  1024  px,  each  pixel  is  square 
with  a  side  length  of  6.45  //m.  Three  60  mm  focal  length  lenses  are  placed  at  the  corners  of 
an  equilateral  triangle  inscribed  in  a  circle  of  radius  of  42.5  mm.  The  CCDs  are  off-centered 
an  additional  0.84  mm  compared  to  the  calculated  values  in  the  alignment  section  so  that 
the  entire  CCD  is  used.  This  increases  the  number  of  pixels  used  by  58%  relative  to  the 
setup  without  the  offset.  The  target  plate  used  for  alignment  was  placed  625  mm  from  the 
camera.  The  2:-range  of  the  test  volume  is  517  mm  to  567  mm  from  the  camera. 

Using  all  the  results  that  were  obtained  before,  the  sensitivities  in  x  and  y  range  from  18.1  to 
19.9px/mm  and  the  sensitivity  in  2:  ranges  from  2.35  to  2.83px/mm,  within  the  test  volume. 
The  ratio  of  the  sensitivities  is  therefore  between  7.0  and  7.7. 
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3.3  Data  Processing  Program 


The  data  processing  program  was  written  to  convert  a  set  of  three  images  into  3D  positions 
and  when  required,  computes  a  velocity  field  from  those  results. 

The  images  are  first  Fourier  transformed  to  remove  noise  inherent  to  the  camera  system. 
Next  a  smoothing  filter  blurs  the  image  with  a  Gaussian  blur.  By  doing  this  the  blurred 
particles  remain  Gaussian  and  their  centroid  does  not  change.  After  the  additional  blurring 
the  noise  is  greatly  reduced.  The  final  step  is  to  scale  the  images,  by  matching  the  contrast 
of  the  three  images. 

After  these  image  manipulations,  the  locations  of  the  peaks  on  each  image  are  searched  for. 
The  positions  of  the  peaks  are  corrected  by  applying  the  calibration  results.  Using  the  three 
CCD  images,  triangles  are  formed  from  corresponding  peaks.  Three  dimensional  particle 
positions  can  then  be  calculated  and  particle  tracking  is  performed  between  two  successive 
frames. 

The  data  processing  program  will  be  explained  in  more  detail  in  a  future  publication  [7]. 


3.4  Error  Analysis 

Without  a  thorough  understanding  of  the  errors  associated  with  3DPTV,  the  results  have 
very  little  meaning.  It  is  therefore  important  to  understand  the  many  different  types  of 
errors.  Because  of  this  a  thorough  error  analysis  was  performed,  both  theoretically  and 
experimentally. 

The  first  of  these  errors  is  the  manufacturing  error.  This  error  takes  into  account  all  possible 
inaccuracies  in  the  lens  positions  as  well  as  inaccuracies  in  the  physical  construction  of  the 
calibration  setup.  These  errors  lead  to  an  error  estimation  which  is  proportional  to  the 
measured  velocities.  The  second  error  is  due  to  the  data  processing  program.  Here,  the  error 
is  constant  for  all  velocities  and  depends  mainly  on  the  depth  position  of  the  particle.  A 
thorough  discussion  of  these  errors  is  out  of  scope  for  this  article  but  will  be  discussed  in  a 
future  paper [8]. 

In  addition  to  the  errors  of  the  3DPTV  system,  other  inaccuracies  have  to  be  considered.  In 
a  supersonic  flow,  changes  of  indices  of  refraction  have  to  be  considered  when  evaluating  the 
3D  position  of  each  particle.  In  addition,  the  finite  size  of  particles  is  important  wherever 
high  velocity  gradients  exist.  These  inaccuracies  will  be  discussed  in  this  section. 


3.4.1  Index  of  Refraction 

The  fluid  inside  and  outside  the  test  section  is  subject  to  change  in  its  refractive  index  since 
the  fluid  density  is  not  constant.  The  refractive  index  of  air  is  roughly  a  linear  function  of 
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its  density 


n  =  l  +  K-^, 

Prei 

where  K  =  0.000292  and  pref  =  1.225kg/m^.  In  the  Ludwieg  Tube  facility,  with  an  initial 
pressure  of  100  kPa,  the  density  could  vary  between  0.15  kg/m^  and  1  kg/m^  and  therefore, 
the  refractive  index  can  vary  between  1.00004  and  1.0002.  The  refractive  index  remains  very 
close  to  1  in  all  cases.  It  is  also  important  to  note,  that  even  if  the  indices  of  refraction  were 
constant  they  would  not  be  known  exactly,  and  hence  would  cause  an  error. 

Index  of  Refraction  uncertainty 

The  indices  of  refraction  of  the  fluid  and  the  glass  are  not  known  exactly;  therefore,  the 
actual  particle  location  can  not  be  calculated  from  the  apparent  particle  location.  Assuming 
we  have  three  indices  of  refraction  (e.  g.,  test-section  gas,  window,  and  outside  air)  and  that 

the  angles  are  small,  the  error,  e,  on  the  CCD  due  to  uncertainties  in  the  indices  of  refraction 
is  given  by 


where 

h  =  ^hi+hi, 

hx  and  hy  being  the  x-  and  y-distances,  respectively,  from  the  lens  center  to  the  image 
location  on  the  CCD.  ni  and  ns  are  the  indices  of  refraction  of  the  fluid  inside  and  outside 
of  the  test  section,  respectively,  while  n2  is  the  index  of  refraction  of  the  glass.  An^  is  the 
uncertainty  in  the  index  of  refraction  n*,  is  the  distance  from  the  camera  to  the  glass, 
is  the  thickness  of  the  glass,  and  di  is  the  distance  from  the  glass  to  the  particle. 

Finite  Density  Gradients 


In  the  test  section  of  supersonic  flow,  significant  finite  density  gradients  will  eixist.  These 
density  gradients  will  cause  a  change  in  the  index  of  refraction  of  the  fluid.  Assuming  that 
the  index  of  refraction  of  all  three  media  are  approximately  equal,  one  can  find  the  error 
in  terms  of  the  index  of  refraction  gradient,  y.  It  is  also  assumed  that  the  particle  is  at 
a  distance  2  from  the  lens  and  that  it  is  inside  the  density  gradient,  that  the  gradient  is 
perpendicular  to  the  light  path,  and  that  the  light  travels  through  the  density  gradient  for 
a  distance,  Lg.  The  error,  e,  on  the  CCD  is 

e  _  dnTg 
dc  y  2z' 


Infinite  Density  Gradients 

Shocks  produce  what  can  be  considered  infinite  density  gradients.  As  in  the  case  of  a  finite 
density  gradient,  these  density  gradients  create  a  change  in  the  index  of  refraction  of  the 
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fluid.  We  will  assume  that  the  index  of  refraction  is  approximately  constant  everywhere,  and 
that  there  is  only  a  change  in  the  index  of  refraction  across  the  shock.  The  angle  between 
the  light  path  and  the  density  gradient  is  a.  The  distance  from  the  particle  to  the  shock  is 
I/s  and  the  distance  from  the  particle  to  the  lens  is  2:.  The  error,  e,  on  the  CCD  is  then  given 
by 


e  ,  An  Z/g 

—  =  tan  (a) - . 

dc  n  z 


3.4.2  Tracking  Error 

A  significant  error  in  the  3DPTV  system  is  not  caused  by  the  camera  itself,  but  rather  the 
ability  of  the  seeding  particles  to  accurately  follow  the  flow.  Because  the' particles  have  mass, 
they  are  unable  to  respond  instantaneously  to  the  flow.  Most  particles  suitable  for  seeding 
supersonic  flow  have  a  speciflc  weight  near  one. 

The  speed  of  a  particle.  Up,  in  Stokes  flow,  across  of  shock  with  speed  Ui  before  the  shock, 
and  speed  U2  behind  the  shock  is  ’ 

—  + 

Up  —  U2  +  {Ui  —  U2)  c  . 

Therefore  the  time  response  across  a  shock  is  characterized  by  a  time  constant,  r,  which  is 
defined  as 

9// 

As  could  be  expected,  the  time  constant  depends  only  on  the  particle  radius,  the  particle 
density,  and  the  viscosity  of  the  fluid.  Given  that  the  viscosity  is  fixed  by  the  choice  of 
flow  conditions,  should  be  minimized.  However,  since  the  amount  of  light  incident  on  the 
particle  is  proportional  to  this  would  also  decreases  the  scattered  light  which  is  used  to 
detect  particles.  The  time  constant  also  decreases  with  decreasing  particle  density. 

Because  of  the  momentum  of  the  particles,  the  particles  will  form  what  can  be  considered 
as  an  apparent  shock  thickness.  This  shock  thickness  can  be  defined  as  the  region  around 
a  shock  where  the  velocity  of  the  particle  differs  from  the  flow  velocity  by  more  than  some 
percentage.  Since  between  the  two  frames  the  particle  moves  some  finite  distance,  3DPTV 
is  not  a  point  measurement  system.  Rather,  the  velocity  measurement  is  an  integration  of 
the  velocity  from  start  to  finish.  Because  of  this,  it  is  not  clear  where  the  velocity  vector 
should  be  placed;  therefore,  the  current  authors  chose  to  place  the  velocity  vector  at  the 
midpoint  of  the  start  and  end  positions.  This  placement  becomes  particularly  important 
when  examining  the  velocity  error  around  a  shock.  For  Stokes  flow,  there  is  an  analytic 
solution  for  the  apparent  shock  thickness  for  a  given  velocity  tolerance.  Figure  9  shows  the 
relationship  between  the  shock  thickness  and  the  velocity  tolerance.  There  are  three  distinct 
cases  that  need  to  be  considered.  The  first  case  is  where  a  particle  starts  upstream  of  the 
shock  and  finishes  downstream  of  the  shock,  but  the  midpoint  between  the  start  and  end 
positions  is  upstream  of  the  shock,  this  case  forms  the  curve  left  of  the  first  dashed  line. 


17 


Apparent  Shock  Thickness 

Figure  9:  Apparent  shock  thickness  due  to  finite  mass  of  particles. 


The  second  case  is  similar  to  the  first  case,  except  that  the  midpoint  between  the  start 
and  end  positions  is  downstream  of  the  shock,  and  this  case  forms  the  curve  between  the 
two  dashed  lines.  The  final  case  is  when  the  start  and  end  positions  of  a  particle  are  both 
downstream  of  the  shock,  and  forms  the  remainder  of  the  curve.  The  theory  supposes  that 
the  particle  follows  the  flow  perfectly  upstream  of  the  shock.  We  can  see  from  the  figure 
that  the  upstream  apparent  shock  thickness  is  small  compared  to  the  downstream  thickness 
when  the  time  duration  between  frames  is  of  the  order  of  the  time  constant.  Because  of 
this,  it  is  reasonable  to  consider  the  downstream  thickness  as  the  total  shock  thickness.  The 
downstream  shock  thickness,  4,  is  given  by 


U2T 


where  f/  —  ^  —  l,t— ^,4  being  the  time  between  frames,  and  e  is  the  velocity  ratio, 
used  to  define  the  shock  thickness  with  U  being  the  velocity  of  the  flow. 

Discretization  Error 


u  ’ 


Part  of  the  error  in  peak  searching  is  due  to  discretization  error.  A  theoretical  approximation 
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for  the  discretization  error,  for  rb  >  1,  shows  that  the  error,  e,  is 

1 


where  A  is  the  amplitude  of  a  particle  image  on  the  CCD.  For  a  given  particle  brightness  the 
product  Arl  is  constant.  Therefore  for  a  constant  particle  brightness  the  error  is  proportional 
to  the  blur  radius.  The  constant  of  proportionality  is  best  found  by  numerical  simulations 
for  a  particular  peak  searching  algorithm. 

CCD  Noise  Error 

A  theoretical  approximation  for  the  error  due  to  noise  on  the  CCD,  e,  can  be  found.  It 
proves  that,  for  rb  >  1,  the  error  is 

E 

(6) 

where  E  is  the  amplitude  of  the  noise.  For  a  constant  particle  brightness  the  error  is  propor¬ 
tional  to  the  blur  radius  squared.  Once  again,  the  constant  of  proportionality  can  be  found 
by  numerical  simulations. 


3.5  Results 

In  this  section,  several  results  will  be  presented.  The  first  of  these  are  theoretical  results 
to  examine  the  highest  accuracy  that  can  be  obtained  using  the  program.  A  moving  plate 
is  then  examined  to  determine  the  errors  one  can  get  on  displacements  and  positions.  A 
slightly  less  serious  but  quite  impressive  result  is  then  presented,  which  shows,  using  a  doll 
face  how  well  a  surface  can  be  represented  in  a  single  measurement.  Finally,  a  supersonic 
flow  result  will  be  presented. 


3.5.1  Synthetic  Images 

An  important  test  of  the  3DPTV  system  was  the  processing  of  synthetic  images.  These 
images  allow  for  a  complete  testing  of  the  program,  including  triangle  formation  and  trian¬ 
gulation  calculations. 

The  domain  was  divided  into  21  planes  perpendicular  to  the  camera  axis.  In  each  of  these 
planes,  400  particles  were  created  and  their  images  were  added  to  form  the  synthetic  images. 
Therefore,  each  synthetic  image  contained  8,400  particle  images.  Each  particle  image  was 
Gaussian  blurred  with  a  10%  blur  radius  of  approximately  3  px.  A  total  of  100  sets  of  images 
were  created  and  then  processed.  In  each  of  the  planes,  the  results  were  averaged  to  evaluate 
the  errors  (*.  e.,  each  symbol  on  Figures  10  and  11  represents  an  average  of  40,000  results). 

As  expected,  the  mean  error  tends  to  zero.  In  each  of  the  planes,  the  mean  error  is  lower 
than  0.009  ^m  and  0.09  ^m  in  the  x-  and  y-directions  and  in  the  ^-direction,  respectively. 


19 


z  position  (in  mm) 


Figure  10:  90%  deviation  error  of  particle  positions  from  synthetic  images 

The  90%  deviation  error  was  also  evaluated  and  is  shown  in  Figure  10.  As  expected,  this 
error  varies  according  to  theory  [8]. 

The  ratio  of  the  90%  deviation  errors  was  also  evaluated,  as  is  shown  in  Figure  11,  and  shows 
a  very  good  agreement  with  the  trend  line  obtained  theoretically  [8]. 

3.5.2  Moving  Plate 

The  moving  plate  test  case  was  performed  using  a  flat  plate  with  an  array  of  white  dots 
on  a  black  background.  There  were  100  dots  per  square  inch.  The  plate  was  placed  at 
z-positions  of  469,  494,  519,  529,  539,  549,  559,  569,  594,  and  619  mm  from  the  camera. 
The  plate  was  then  moved  in  three  equal  increments  of  2.54  mm  in  the  y-direction  for  each 
of  the  ^-locations,  this  resulted  in  a  total  of  6  displacement  vectors  for  each  2:-position  (3 
displacements  of  2.54  mm,  2  displacements  of  5.08  mm  and  1  displacement  of  7.62  mm). 
Similarly  the  plate  was  moved  in  three  equal  increments  of  2.54  mm  in  the  ^-direction  for 
each  z-position. 

First,  the  mean  error  was  evaluated.  This  error  is  proportional  to  the  displacement.  Figure  12 
shows  the  relative  mean  error  obtained  for  2.54  mm  and  5.08  mm  displacements,  each  data 
point  is  an  average  of  3  and  2  values,  respectively.  In  this  plot,  the  7.62  mm  results  are  not 
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Figure  12;  Mean  plate  displacement  error.  White,  grey,  and  black  symbols  refer  to  the  errors 
in  the  x-,  y-,  and  ^-directions  respectively.  The  2.54  mm  and  the  5.08  mm  cases  are  the 
average  of  3  and  2  data  points,  respectively. 

shown  since  they  correspond  to  the  average  of  the  three  2.54  mm  results. 

Even  though  the  mean  errors  are  very  small,  all  within  2%,  they  remain  in  general  far 
above  the  results  expected  theoretically  [8].  One  possible  reason  for  the  inaccuracies  of  the 
y-measurement  in  the  .^-displacement  case,  is  that  the  whole  camera  was  moved  instead  of 
the  plate.  Due  to  the  weight  of  the  camera  and  the  lack  of  robustness  of  the  ym  stages,  it 
might  have  changed  its  pitch  angle  very  slightly  when  moved.  An  inaccuracy  of  20  y,m  per 
mm  in  the  y-measurement  can  be  obtained  by  having  a  change  in  pitch  angle  of  less  than 
5  X  10~®rad  per  mm  (less  than  0.003°  per  mm).  Furthermore,  the  resolution  of  the  gauges 
that  measured  the  displacements  is  25.4  ym,  it  is  possible  that  part  of  the  error  comes  from 
inaccuracies  in  the  displacement  of  the  plate. 

The  90%  deviation  error  of  the  displacements  was  also  evaluated.  The  results  are  plotted  in 
Figure  13.  Each  symbol  is  the  average  of  the  90%  deviation  errors  from  three  and  two  data 
sets  in  the  2.54  mm  and  the  5.08  mm  cases,  respectively. 

As  expected  from  theory,  the  error  can  be  separated  into  two  parts.  The  first  of  these  is 
proportional  to  the  displacement  and  the  second  is  constant  for  all  displacements.  The  total 
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Figure  13:  90%  deviation  errors  of  plate  displacement.  White,  grey,  and  black  symbols  refer 
to  the  errors  in  the  x-,  y-,  and  2;-directions  respectively.  The  2.54  mm  and  the  5.08  mm  cases 
are  the  average  of  3  and  2  data  points,  respectively. 
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Figure  14:  90%  deviation  errors  of  plate  displacements.  White,  grey,  and  black  symbols  refer 
to  the  errors  in  the  rr-,  y-,  and  ^-directions,  respectively. 


90%  deviation  error,  6D^  can  therefore  be  written  as 

SD=^ClD^  +  CI.  (7) 

Since,  for  most  flow  measurements,  the  displacement  of  the  particles  will  be  smaller  than 
2.54  mm,  it  is  interesting  to  analyze  the  results  for  this  particular  displacement,  as  this 
represents  the  worst  case  scenario.  Figure  14  therefore  gives  a  higher  bound  for  the  90% 
deviation  error  one  can  expect  with  the  3DPTV  setup. 

According  to  these  results,  the  90%  deviation  error  in  the  test  volume  remains  under  5.2yLim 
and  under  63//m  in  the  x-,  ^/-directions  and  in  the  2;-direction,  respectively. 

Both  the  theoretical  values  and  the  experimental  data  of  the  90%  displacement  error  are 
shown  in  Figure  15.  The  discrepancy  in  the  results  for  low  2^-position  comes  from  the  fact 
that  the  dots  on  the  plate  appear  very  blurred  and  that  the  amplitude  of  the  peaks  becomes 
very  low.  This  effect  is  presented  in  Equations  5  and  6. 

The  varying  part  of  the  90%  deviation  error  is  presented  in  Figure  16.  Although  these  errors 
remain  quite  small  (lp,m  per  mm  in  a;  or  y  and  10  nm  per  mm  in  z),  they  are  much  higher 
than  what  was  expected.  A  possible  reason  for  that  is  again  the  lack  of  robustness  of  the 
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Figure  15:  Fixed  part  of  90%  deviation  errors  of  plate  displacements.  White,  grey,  and  black 
symbols  refer  to  the  errors  in  the  x-,  y-,  and  2:-directions,  respectively. 
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Figure  16:  Variable  part  of  90%  deviation  errors  of  plate  displacements.  White,  grey,  and 
black  symbols  refer  to  the  errors  in  the  x-,  y-,  and  ^-directions,  respectively. 

translation  stages.  In  fact,  if  during  the  movement  of  the  plate  in  the  y-direction,  its  pitch 
angle  changed  by  only  4.5  x  10“'^  rad  per  mm  (0.026°  per  mm),  the  90%  deviation  in  ^  would 
reach  10  /rm.  A  small  roll  angle  of  4.5  x  10~®  rad  per  mm  (0.0026°  per  mm)  is  sufficient  to 
explain  the  inaccuracies  in  x  and  y. 

^From  all  the  results  presented  here,  it  is  possible  to  obtain  an  upper  bound  for  the  velocity 
errors.  The  mean  errors  will  be  approximately 
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whereas  an  estimation  of  90%  deviation  errors  of  the  velocities  is 


5U,= 

6Uy  = 

5U,= 


1  c:  V  in-4 

'751/2  +  12C/2  +  gij2  + . 0-4 

/  1 ^  V  1  n-4 

'12//2  +  75C/2  +  9172  +  — . /  X  10"^ 


I87U^  +  87t/2  +  0.85t/2  +  x  10“^ 

4 


(11) 

(12) 

(13) 


where  the  velocities  and  the  time  between  frame,  ts,  are  in  m/s  and  in  s  ,  respectively. 

If  in  fact,  most  of  the  disagreement  with  theory  is  due  to  the  inaccuracy  of  the  plate  place¬ 
ment,  these  errors  could  be  much  smaller.  Since  the  purpose  of  the  present  work  is  not  to 
find  the  exact  accuracy  of  the  system  and  that,  due  to  inertia,  the  actual  error  will  depend 
most  heavily  on  the  particles  ability  to  track  the  flow.  Equations  8  through  13  are  only  an 
estimate  of  the  velocimetry  errors. 


3.5.3  3D  Solid  Imaging 

An  additional  capability  of  the  3DPTV  system  is  the  ability  to  do  3D  solid  imaging.  To 
illustrate  this  capability  the  system  was  used  to  image  the  face  of  a  doll.  This  was  done  by 
painting  a  doll  face  black,  and  putting  a  large  number  of  white  dots  on  the  face.  While  it  is 
hard  to  make  analytic  conclusions  from  the  results.  Figure  17  shows  excellent  resolution  of 
the  face.  Some  of  the  roughness  on  the  image  is  possibly  due  to  the  fact  that  some  of  the 
white  dots  were  large  and  non-circular.  Note  that  these  results  were  obtained  using  an  old 
version  camera,  which  did  not  take  advantage  of  offsets.  This  camera  was  able  to  investigate 
a  10  X  10  X  10  cm^  volume  but  with  a  lower  accuracy. 


3.5.4  Sonic  Jet  Injection 

The  ultimate  goal  of  the  current  authors  is  to  use  the  3DPTV  system  to  image  supersonic 
flows.  Particle  tracking  measurements  in  supersonic  flows  are  particularly  challenging  and 
present  many  problems  not  seen  in  slower  flows.  The  primary  problem  is  that  in  order  for 
the  particles  to  respond  quickly  to  the  flow  they  must  be  very  small. 

Experiments  using  a  sonic  jet  injector  with  1  /rm  diameter  olive  oil  droplets  were  conducted. 
Even  with  a  200  mJ  Nd:YAG  laser  these  particles  were  very  dim  and  detection  was  very 
difficult.  Despite  the  low  signal  to  noise  ratio  the  3DPTV  system  was  able  to  produce  an 
accurate  velocity  field.  Figure  18  shows  an  instantaneous  velocity  field  containing  224  vectors, 
which  is  approximately  35  vectors  per  cubic  centimeter.  The  particles  were  accelerated  to 
sonic  speeds  through  a  converging  nozzle  with  an  outlet  diameter  of  7  mm.  The  region 
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Figure  18:  Velocity  vectors  of  a  sonic  jet  injector  with  1  /rm  olive  oil  seeding  particles. 

downstream  of  the  injector  nozzle  was  kept  in  vacuum;  therefore,  the  particles  continued  to 
accelerate  to  supersonic  speeds  after  exiting  the  injector. 


3.6  Conclusions 

The  3DPTV  system  has  been  shown  to  have  excellent  accuracy.  For  velocities  parallel  to  the 
camera  face  the  accuracy  is  within  0.25%  and  for  velocities  perpendicular  to  the  camera  face 
the  accuracy  is  within  2.5%.  These  errors  are  well  understood  and  agree  with  the  theoretical 
estimates.  Solid  surface  imaging  of  a  doll  face  provided  impressive,  albeit  non-quantitative 
results.  Most  importantly,  velocity  measurements  of  a  sonic  jet  were  made.  These  results 
are  promising,  although  the  vector  density  is  lower  than  desired.  Future  work  has  begun 
on  using  larger  and  lighter  particles,  which  will  significantly  improve  the  luminosity  of  the 


particles.  With  brighter  particles  the  results  should  begin  to  approach  that  of  the  synthetic 
images  and  flow  flelds  with  over  5,000  vectors  should  be  obtainable. 


4  General  Conclusions 


The  supersonic  Ludwieg  Tube  at  GALCIT  was  completed  and  taken  into  operation.  The 
injection  Ludwieg  Tube  was  constructed  and  taken  into  operation,  thus  completing  a  com¬ 
pact,  inexpensive  supersonic  injection  facility.  The  implementation  of  three-component, 
three-dimensional  particle  image  velocimetry  in  supersonic  air  flows  from  the  instrument 
existing  for  low  speed  water  flows  at  the  time  of  inception  of  the  program  required  a  much 
greater  development  effort  than  had  been  anticipated  at  the  outset.  Major  redesign  of  the 
hardware  and  new  software  development  as  well  as  experimentation  with  particle-genrating 
systems  produced  an  instrument  that  is  now  able  to  fulfill  the  tasks  demanded  of  it  for  the 
project  with  good  accuracy  and  spatial  resolution. 
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V 


APPENDIX 

Original  Statement  of  Work 
Year  1.: 


1.  Design  and  construct  the  injection  system.  This  consists  of  the  injection  Ludwieg 
tube  (ILT)  leading  to  the  injection  nozzle  in  the  test  section  via  a  quick-acting  valve. 
Design  and  construct  the  injection  wall,  for  generating  a  new,  laminar  boundary  layer 
upstream  of  the  injection  point. 

2.  Shake  down  the  ILT  and  its  synchronized  operation  with  the  LT.  Test  out  parameter 
range  limits  using  schlieren  visualization. 

3.  Design  and  construct  the  DDPIV  system,  including  suitable  optical  access  to  the  test 
section,  purchase  and  construction  of  triple  camera,  purchase  and  assembly  of  comput¬ 
ers  and  frame  grabbers,  installation  of  software. 

4.  Test  and  calibrate  the  complete  system,  shake  down  and  find  limits  of  spatial  resolution 
and  flow  parameter  ranges. 


Years  2  and  3: 


1.  In  consultation  with  the  computational  fluid  dynamics  community,  define  the  test  cases 
of  particular  interest  and  concentrate  initial  efforts  on  those,  in  order  to  establish  a 
rapport  with  the  computational  groups. 

2.  Explore  the  parameter  space  (Reynolds  number,  momentum  ratio,  oncoming  boundary 
layer  state,  injection  Mach  number,  density  ratio)  in  order  to  discover  where  qualitative 
changes  occur.  Delineate  the  boundaries  on  a  map  of  parameter  space. 

3.  Build  database  of  detailed  flow  field  measurements  in  specific  regions  of  the  parameter 
space. 

4.  Document  results  in  easily  accessible  electronic  form  and  make  them  openly  available 
to  researchers  approved  by  the  AFOSR. 

5.  Document  results  in  publications 
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